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EFFECTS OF MAGNETICALLY INDUCED EDDY-CURRENT TORQUES 


ON SPIN MOTIONS OF AN EARTH SATELLITE 1 
By G- Louis Smith 


SUMMARY 


One of the sources of torques on near-earth satellites is the interaction 
of the earth 1 s magnetic field with eddy currents induced in the electrically 
conducting parts of a spinning satellite. Initially, it is assumed that the 
geomagnetic field is that of a space-fixed dipole. The equation for the average 
torque, due to this effect, acting on a satellite of known properties is derived 
and presented. This equation shows the variation of the torque with the orbital 
parameters. The time history of the spin vector as influenced by this torque 
has been investigated. For the general case, the spin vector can be resolved 
into three orthogonal components which are damped exponentially at three differ- 
ent rates. 

The analysis is extended to include the effect of the tilt of the geomag- 
netic dipole axis with respect to the earth* s spin axis. The analysis and equa- 
tions which include this effect closely parallel the preceding analysis and 
equations . 

The generalized parameters, which can be used in analyzing the motion due 
to eddy currents of any spinning satellite in a near-circular orbit, are pre- 
sented in graphical form. 


INTRODUCTION 


One of the sources of torques on near-earth satellites is the interaction 
of the earth* s magnetic field with eddy currents induced in the electrically 
conducting parts of a spinning satellite. This effect is considered in refer- 
ences 1 to 6. The problem consists of two parts: first, the calculation or 

measurement of the eddy-current torques on the satellite due to the presence of 
a given magnetic field, and second, the calculation of the torques and their 
effect on the angular motions of the satellite while in orbit around the earth. 
Studies of the first part of the problem are presented in the literature, but 
the latter part of the problem is treated only by estimates or else for very 
special cases. Sanduleak (ref. 7 ) sets up the equations for the spin-damping 
torque, in which the earth* s magnetic field and orientation of the spin vector 
are considered as variables. Sanduleak programed these equations on an elec- 
tronic data processing system for the torque computation, and numerical integra- 
tion of the equations of motion yields the spin-rate history. However, only the 
iThis "Corrected Copy* 1 of NASA TN D -2198 includes corrections to pages 7, 
15, and l6 and replacement of figures 2, 3, 7 * 8, and 9 * An appendix is 

added to clarify the derivation of equation (l*0, which was the source of the 
errors . 



damping component of torque is considered. It is shown by Vinti in reference k 
that there is also a torque component tending to precess the spin vector. 

\ 

The purpose of this paper is to derive and present expressions for the 
average torques due to induced eddy currents acting on a satellite and to study 
the resulting time history of the spin motion. Examples of application of this 
study are the improved estimation of torques acting on a spinning satellite 
which must be oriented or for which the spin must be maintained within limits, 
and the calculation of the damping of satellite spin for comparison with 
observed variations. 
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A Cartesian coordinate system is defined; with the Z-axis along the dipole 
axis and the X-axis passing through the ascending node of the intersection of 
the orbit plane with the geomagnetic equatorial plane. The Y-axis then completes 
the coordinate system; as shown in figure 1. The orbit orientation is described 
with respect to this system by the inclination to the geomagnetic equator i 
and the argument of the perigee q. The position of the satellite is given by 
the orbit angle 0 or by the longitude 0 and colatitude i|r within the coor- 
dinate system. 

The torque due to eddy currents acting on a satellite is shown by Vinti 
(ref. 4) to be 



(i) 


if 



where K is a constant for the satellite spinning about a given body axis. 
(Small variations of K due to such effects as change of conductivity with 
satellite temperature are neglected in this paper.) Equation (l) is rewritten 
as 


( 2 ) 


The time required for the spin to decrease to one-half its original value due 
to eddy-current torques is in general on the order of a hundred days. (See, 
for example, refs. 1 to 3.) Hence, during a day, the change in spin rate is on 
the order of 1 percent, and during a single orbit, approximately 2 hours, the 
spin vector is very nearly constant. If the spin vector is assumed constant, 
the torque may be integrated with respect to time around one orbit, and an 
average torque is thus calculated: 


Mav = i / M dt = | [ 1 *] <(>} (3) 

° o 



where 



It is next necessary to evaluate the elements of 



Evaluation of Matrix Elements 

The geomagnetic field is assumed to be a dipole field and the tilt of the 
axis of the dipole with respect to the spin axis of the earth is neglected; 
therefore, the field is not a function of longitude and time. The magnetic 
field is thereby a function of the position of the satellite, which is given by 
the classical orbit relations. 

The spherical components of a dipole field of strength p are 
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By orbit mechanics 
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1+e cos 0 


(8) 


The angles >|r and 0 are given by 


and 


cos i|r = sin i sin(0 + t)) 


tan 0 = cos i tan(0 + T)) 


(9) 

( 10 ) 


The variable of integration in equation (3) is changed to 0j thus. 
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and 
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Equations ( 8 ), ( 9 ), and (10) are used, with the manipulation presented in 
the appendix, in equation ( 6 ) to express the components of if as functions of 
0 and the orbit parameters p, e, i, and tj. These expressions are substi- 
tuted into equation (12) and the integration performed. It is found that terms 

involving e and vanish. If the analysis is restricted to near-circular 

orbits, the terms in c 2 and e + can be neglected, and is found to be 
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Figure 2 shows the elements of [a( i )] as a function of inclination i. 
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The orbital period is replaced by 


and the e 2 term is 


neglectedj thus, 
becomes 
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equation ( 3 ) for the average torque due to eddy currents 


- - ^[A(D] H 
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(15) 


This torque can be added to other torques , for example, those due to grav- 
itation gradients, in order to study the motion of a satellite. In the present 
paper, however, the only external torque considered is that due to eddy currents. 

It is well-known that for spinning satellites, because of internal dissipa- 
tion of kinetic energy of rotation, the only stable axis of rotation is the axis 
of maximum momentum. After the transition to this spin motion is complete, the 

angular momentum of the satellite is simply I mx , and the equation of the 

spin motion of the satellite becomes 


^max ~ ~ (l^) 


The solution form 


^y - ^<&y ^ max ^ 

is substituted into equation (l6), from which 

[ MD ] {^} 0 = ^(“} 0 

or 

([A(D] - a[i]){o>> o - 0 ( 17 ) 

which is simply a matrix eigenvalue problem. Equation (17) has nontrivial 
solutions only if 


8 



a 


0 


0 


0 


(18) 


det 10**0 - "01 


e - >- 


7 - A 


= 0 


It is noted that the x-component is uncoupled from the y- and z-components, 
that is, all terms in the first column and first row vanish, except a - A, and 
the problem is simplified considerably. The three solutions of equation (l8) 
are 
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A 2,3 


P + 7 - 
2 + 



P + 7 
2 + 



(19) 


The eigenvalues are plotted in figure 3 as a function of inclination. 

The components of the corresponding unit eigenvectors are the minors of 

the determinant 11**3 - *QQ| obtained by moving across a row and normal- 
izing and are 
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Tlie matrix ([AUO - X[l]) is symmetric; hence the eigenvectors 


are 


orthogonal and the calculations are simplified considerably. Since e^_ coin- 


cides with the X-axis, e 2 and e^ are obtained from 


j and k by rotating 
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as shown in figure 4 . Thus 
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Comparison of equations (20) and (21) shows that 


( 21 ) 


tan X = 
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This angle of rotation X is shown as a function of orbital inclination in 
figure 5- Finally, the spin vector may be written as 
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where o^q, ay^ 0 , a nd are components of 

coordinate system. 


( 25 ) 


in the Cartesian 


Thus far, it has been assumed, for simplicity, that the geomagnetic field 
is a dipole with the axis along the earth’s spin axis. Consideration is now 
given to the effects of the tilt of the dipole axis with respect to the earth’s 
spin axis. 


Rotating Earth With Tilted Dipole Field 

The geomagnetic field of the earth is approximated by a dipole field, the 
dipole axis being tilted 17° with respect to the earth’s spin axis. (See 
ref. 70 

In order to include the effect of the tilt of the magnetic dipole field 
with respect to the earth’s spin axis, it is necessary to define three coordi- 
nate systems. (See fig. 6.) The space-fixed or inertial system Xj,Yj,Zj is 

oriented with the Zj-axis along the earth’s spin axis, and the Xj-axis along 

the intersection of the satellite orbit plane and the equator so that the 
Xj-axis passes through the ascending node. The earth-fixed system X;e,Y;e,Z;e 
is oriented with the Z E -axis along the earth’s spin axis and the Xg-axis passing 

through the intersection of the equator with the geomagnetic equator. Finally, 
the geomagnetic system X^,Y^,Z^[ is defined with the Z^-axis along the axis 

of the geomagnetic dipole and the X^-axis passing through the ascending node of 

the orbit referenced to the geomagnetic equator. In addition, several angles 
are defined: ij denotes the orbital inclination with respect to the earth’s 

equator and i^ the orbital inclination with respect to the geomagnetic equa- 
tor, 1 jl is the angle between the Xj- and x E -axes, v the angle between the 
X^- and Xg-axes, and finally £ is the tilt of the dipole field with respect 
to the earth’s spin axis, that is, the angle between Z E and Z^. It is seen 

in figure 6 that these angles are all related, in the spherical triangle ABC. 

The purpose of each coordinate system is as follows: 

The inertial coordinate system provides a necessary fixed reference, as the 
other coordinate systems rotate with the earth. The analysis of the previous 
section was referenced to the geomagnetic coordinate system. The earth-fixed 
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system serves as an intermediate coordinate system by which Euler angle rota- 
tions are made in transforming from the inertial to the geomagnetic coordinate 
system, and is further necessitated by the requirement that dp/dt, which is 
the spin rate of earth, is constant. 

As is seen from figure 6, the X E ,Y E ,Z E system is obtained by rotating 
the Xj,Yj,Zj system about the Z-axis through the angle p, or 
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The Xj 4 ,Yjv[,Zf 4 system is next obtained by rotating the X E ,Y E ,Z E system, 
first about the Xg-axis by the tilt angle f; , and then about the Zp[-axi s by 
the angle v : 
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Equations (26) and (27) yield 


X M " 




(28) 


The components of the vector in the inertial coordinate system are given in 
terms of the components in the geomagnetic coordinate system by simply 
reversing the rotations and their order: 

X I = HH[ T -v]x M (29) 

Equation (15) was derived on the assumption of a dipole field constant in 
time. As the earth rotates, the geomagnetic axes also move and therefore equa- 
tion (15) does not apply. If the variation of the dipole field is slow, 


12 



equation (15) may used as an approximation to the torque. During a single 
orbital period of the satellite, the change in the geomagnetic field at a point 
is considered to be small. 

The following approximations are now made : 

(1) Equation (15) is applicable 

(2) The summation of angular impulses over 1 day may be approximated by 

integration 

(3) The satellite spin vector does not change significantly in 1 day 

In order to use equation (l5)j it is necessary only to reference the 
vectors to the geomagnetic coordinate system: 


T v][ T JH(Mav>i = -^[ A ( i M)][ T v][ T J[ T n] {“} 


(30) 


from which 
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where 


is given by 


lli.*)] 11 
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and is seen to be simply the similarity transform of J~A ^i]si) J i* rom the geomag- 
netic coordinate system to the inertial coordinate system. The average daily 
torque is thus written as 


The matrix 


{Mav} j ^ = "p"[ B ( il )] 

[ B (il)] which now takes the place of [A(i„)] is 

[^l)] * h /[K 1 !.*! 


(33) 


(34) 

13 


1 



The integration indicated in equation (3*0 was accomplished numerically. The 


matrix 




being a similarity transform of 


is also symmetric; hence its integral, 


[ B W] 
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, which is symmetric. 


, is symmetric. The b^ 2 an & 


b^ (hence h2i and b^, by symmetry) terms are found to be only a fraction 

of a percent of the larger terms in the matrix , and may be considered to be 
negligible. The remaining predominant matrix elements bp^, ^22 j b^j, an< ^ 

b23 (=^32) are plotted as functions of ij, the inclination of the orbit plane 

to the equator, in figure 7- 


Equation (33) can be used to calculate the daily average torque due to 
eddy currents on a spinning satellite. The time history of the effect of this 
torque can be determined by the same analysis as used in the previous section; 


equations (15) to (25) are immediately applicable when the elements of 


are substituted for the elements of 


[*(*>]■ 


— * — 


The eigenvalues of 


B(il) 

- — 


V l> 


A2, and Aj are shown in figure 8 as a function of ij, and the angle Xj 


between the eigenvectors of 



and the space-fixed coordinate axes is 


shown in figure 9 also as a function of ij. 


Numerical Example 

In order to illustrate the application of the theory developed in this 
paper, the spin decay of a typical satellite is calculated. For the calcula- 
tions, the following numbers are used: 


and 


Argument of perigee = 142.3° 
I = 2.0 X 10T gm-cm^ 

K = ,50 dyne-cm-sec/gauss^ 

p = 1.15R e 
ij = 52° 

s = 0.3151 


which is normalized with respect to Rgl The spin vector is considered to 
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be tangent initially to the orbit at the perigee, but in opposite direction 
from the velocity vector. 

The initial spin vector, expressed in the space-fixed-coordinate system, 
is thus 
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It is seen from figure 9 that for ij - 52°, Xj = 33°; thus, by equation (25) 
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From equation (24), 
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The eigenvalues for ij = 52° are shown by figure 8 to be 

A L = 1.20 

A2 = 0.92 
Aj = 1.69 

By using equation (21) for the eigenvectors, the equation for the spin vector 
(eq. (23)) thus becomes 
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CONCLUDING REMARKS 


An analysis has been made to determine the average, or net, torques due to 
magnetically induced eddy currents, which act on a spinning near-earth satellite. 
Initially, the earth’s magnetic field is assumed to be that of a space-fixed 
dipole. An equation for the torque vector, which is a linear function of the 
satellite spin vector, is developed, including the orbital parameters. It is 
found that the first- and third-order terms in eccentricity vanish identically; 
second- and fourth-order terms are neglected. The effect of orbital inclination 
is given by a matrix function, since in general the torque vector is not par- 
allel to the spin vector. 

The time history of the spin vector subject to this torque has been inves- 
tigated. It is found that, for the general case, the spin vector can be 
resolved into three orthogonal components which are damped exponentially at 
three different rates. 

The effect of the tilt of the geomagnetic dipole axis with respect to the 
earth’s spin axis has been considered for both the torque and the resulting 
spin history. The analysis and resulting equations closely parallel those of 
the foregoing analysis which did not consider this effect. The two analyses, 
with and without consideration of the dipole tilt, yield qualitatively similar 
results, although they differ slightly quantitatively, as is to be expected. 

The generalized parameters which can be used in analyzing the motion, due 
to eddy current torques for any spinning satellite in a near-circular orbit are 
presented in graphical form. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., December 10, 1963 . 
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APPENDIX 


DERIVATION OF [A(i)] MATRIX ELEMENTS 


In order to derive the elements of the {a( i J] matrix as expressed in 
equation (l4), equations (6) to (13) are used. In figure 1, note the right 
spherical triangle with sides 0 + t], 0, and ~ - \|r, from which the spherical 

trigonometric relation 
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is obtained. From equation (10), 
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Direct substitution of equations (Al), (A2), and ( 9 ) into equation (6) gives 
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H z 


— (l + e cos 
P 5 



sin^i sin 2 (0 + r|) 
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When equations (A 3 ) are used to expand the integrand of equation (12), a 
fourth-order polynomial in e is obtained. The terms involving e and e5 
have only odd powers of sin 0 and cos 0 and therefore vanish identically 
when integrated from 0=0 to 2jt. The terms involving and do not 

vanish, and the error introduced by neglecting terms in eccentricity is of the 

order of It is noted that the e ^ terms contain sin 2 t) and cos 2 tj« 

If only terms of zeroeth order in e are retained, the integrals of 
interest are 
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where 0 * = 0 + tj. These integrals, with 2 jt factored out in accordance 

]? 

with equation (13), thus give the elements of (^-(1^ : 
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Figure 1.- Coordinate systems and symbols used for analysis of noninclined dipole. XY plane 
is earth equatorial plane, dipole axis is along Z-axis. 
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Figure 3.- Eigenvalues of [Mi)] matrix as a function of orbital inclination 
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Figure 5*- Angle X between eigenvectors and coordinate axes as function of orbital 

inclination i. 
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Figure 6.- Coordinate systems used for calculation of torques due to tilted dipole. 



Figure 7.- Elements of [B(ijY"| matrix as function ij 















"The aeronautical and space activities of the United States shall be 
conducted so as to co?itribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof 

— National Aeronautics and Space Act of 1958 
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